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DEFINITION  OF  TERMS 


Disturbance  Amplification  Ratio 
Skin  Friction  Coefficient 

The  Amplification  Ratio  (=  A/Aq)  according  to  linear  stability  theory 

Disturbance  Frequency  -  Hz 

Boundary  Layer  Shape  Factor 

Reference  Length 

Reynolds  Number  =  U^o  L/^oo 

Displacement  Thickness  Reynolds  Number  =  Uoo^*/i'„ 

Temperature 
Free  Stream  Velocity 

Local  Relative  Growth  Rate  of  a  Disturbance 

Displacement  Thickness 

Free  Stream  Kinematic  Viscosity 

Nondimensional  Disturbance  Frequency  =  27rfi/oo/Uoo^ 
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INTRODUCTION 


This  report  describes  the  implementation  and  evaluation  of  the  Fortran  IV  computer 
program  DABL  (Disturbance  Amplification  in  Boundary  Layers  —  developed  at  the  David 
Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC)  by  vonKerczek  and 
Groves,  1978)  which  calculates  the  growth  of  small  disturbances  in  laminar  boundary  layers. 
DABL  calculates  the  amplification  ratios,  A/Aq,  for  heated  or  unheated,  two-dimensional  or 
axisymmetric,  incompressible  laminar  boundary  layers.  The  boundary  layer  is  assumed 
to  be  steady  and  has  constant  density  with  the  fluid  viscosity  varying  with  temperature.  The 
DABL  code  requires  only  the  body  coordinates,  surface  temperature  distribution  and  body 
Reynolds  number  to  compute  the  spatial  amplification  ratios  as  a  function  of  arc  length  and 
nondimensional  disturbance  frequency  using  parallel,  viscous  flow,  linear  stability  theory. 
The  resulting  amplification  ratio  distribution  can  be  used  to  locate  boundary  layer  transition 
via  the  e^  method. 

This  report  makes  no  attempt  to  discuss  or  evaluate  the  theoretical  and  numerical 
methods  used  by  DABL  in  computing  the  boundary  layer  properties  and  in  performing  the 
boundary  layer  stability  analysis.  vonKerczek  and  Groves  (1978)  discussed  these  topics  in 
detail.  However,  this  report  does  present  several  examples  used  to  evaluate  DABL  as  modi¬ 
fied  by  the  author  and  adapted  to  the  NOSC  UNIVAC  1 1 10  Computer. 

Upon  receipt  of  the  DABL  Program  from  DTNSRDC,  the  computer  code  was 
altered  (to  be  compatible  with  UNIVAC  1110)  and  compiled  in  ASCII  Fortran  utilizing  a 
segmented  execution  mode  to  conform  to  the  65K  storage  limit  of  the  NOSC  system.  Prior 
to  any  further  changes,  DABL  was  checked  for  proper  execution  using  the  examples  of 
vonKerczek  and  Groves  (1978). 

After  the  DABL  code  was  operating  properly,  the  following  additions/changes  were 
incorporated  into  the  code: 

1.  The  program  output  was  expanded  to  include  both  the  original  and  normalized 
body  coordinates  and,  in  the  case  of  heated  bodies,  the  surface  temperature 
distribution. 

2.  The  option  to  read  and  write  the  boundary  layer  data  to  mass  storage  was  added. 

3.  The  fluid  property  definitions  were  expanded  to  include  air,  as  well  as  water, 
boundary  layers. 

4.  Options  were  included  to  allow  the  stability  calculations  to  proceed  regardless  of 
the  calculated  amplification  ratio  value  and  to  begin  the  stability  automatic  start 
routine  at  a  prescribed  body  location. 

5.  The  step  function  surface  temperature  distribution  was  changed  to  allow  the 
distribution  to  end  prior  to  the  end  of  the  body. 

The  remainder  of  the  report  is  divided  into  several  sections.  First,  DABL  was  used  to 
compute  the  neutral  stability  points  for  an  unheated  and  heated  flat  plate  and  the  amplifica¬ 
tion  ratios  at  several  disturbance  frequencies  for  the  unheated  plate.  These  results  are  com¬ 
pared  with  theory  and  experiment.  Next,  a  comparison  is  made  between  the  values  of  skin 
friction,  shape  factor,  growth  rates  and  amplification  ratios  computed  by  DABL  and  TAPS 
(Transition  Analysis  Program  System)  for  both  a  heated  and  unheated  9: 1  ellipsoid. 

Thirdly,  a  complete  description  of  the  input/output  variables  and  program  input  scheme  of 
DABL  is  presented.  Finally,  abbreviated  listings  of  the  computer  output  for  the  heated  flat 
plate  and  the  heated  9: 1  ellipsoid  are  included  as  examples. 
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Although  many  of  the  descriptions  found  in  the  original  DABL  report  (vonKerczek 
and  Groves,  1 978)  are  repeated  here,  the  input  format  for  the  NOSC  version  of  DABL  has 
been  changed  completely.  Anyone  desiring  to  use  this  program  should  contact  the  author 
and  a  listing,  card  deck  or  tape  and  check  solutions  will  be  provided. 


TWO-DIMENSIONAL  GEOMETRY 

To  evaluate  the  two-dimensional  option  of  the  DABL  code,  a  flat  plate  geometry  was 
chosen.  The  test  plate  is  0.25  inches  thick  and  has  a  sharp  leading  edge  with  a  1.8  degree 
half  angle.  There  is  good  agreement  with  the  potential  flow  velocities  computed  by  DABL 
and  the  Douglas  Two-Dimensional  Cascade  Potential  Flow  Program  (Giesing,  1 964).  The 
plate  has  a  nearly  zero  pressure  gradient  at  a  zero  angle  of  attack.  Overall  plate  lengths  of 
one  and  1 2  feet  were  used  in  the  analysis.  The  coordinates  of  the  one-foot  plate  and  the  po¬ 
tential  flow  velocities  are  given  in  Appendix  A. 

The  neutral  stability  curves  for  both  an  unheated  plate  and  a  heated  plate  with  a 
constant  surface  overheat  of  2.78°C  were  calculated  using  the  automatic  starting  option  of 
DABL.  For  each  plate  length  and  Reynolds  number,  the  neutral  stability  points 

(where  -  a*  =  0)  were  located  and  the  corresponding  displacement  thickness  Reynolds  num¬ 
ber  and  nondimensional  frequency  were  tabulated  and  plotted. 

Figure  1  compares  the  neutral  stability  points  computed  by  DABL  (also  found  in 
Table  1)  with  the  experimental  results  of  Schubauer  and  Skramstad  (1948)  and  the  theoreti¬ 
cal  results  presented  by  Ross  et  al  (1970).  As  the  RL/ft  was  increased,  DABL’s  automatic 
starting  method  calculated  nondimensional  frequencies  a>f  which  resulted  in  neutral  points 
falling  within  the  neutral  curve  of  Ross  et  al.  However,  if  for  a  particular  RL/ft,  the  starting 
cjf  is  decreased,  the  DABL  calculated  neutral  points  fall  closer  to  theoretical  and  experimental 
values. 

In  Figure  2,  a  comparison  is  made  between  the  amplification  ratios  A/Aq  computed 
by  DABL  and  the  theoretical  curves  of  Jordinson  (1970)  and  the  experimental  data  presented 
by  Ross  et  al  ( 1 970).  Agreement  is  good  within  the  common  regions  of  comparison.  The 
additional  amplification  ratios  for  Rl/A  ~  2.0  X  10^  and  5.0  X  10^  have  been  included  as 
further  demonstration  of  DABL’s  results;  however,  no  comparisons  are  made. 

Finally,  the  neutral  stability  points  for  a  heated  plate  (constant  AT  =  2.78C)  were  cal¬ 
culated  using  DABL  (also  found  in  Table  2)  and  are  compared  to  the  experimental  results  of 
Strazisar  ( 1 976)  in  Figure  3.  As  in  the  unheated  case,  DABL’s  automatic  starting  mode  com¬ 
putes  frequencies  cuf  which  result  in  neutral  points  falling  within  the  neutral  zone  defined  by 
Strazisar.  In  general,  DABL  produces  conservative  estimates  for  the  onset  of  disturbance 
amplification. 
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Figure  2.  Comparison  of  the  disturbance  amplification  ratios  of  a  flat  plate. 
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Figure  3.  Comparison  of  the  neutral  stability  points  of  a  flat  plate 
with  2.78°C  surface  overheat. 


AXISYMMETRIC  GEOMETRY 

A  9: 1  ellipsoid  was  selected  to  provide  a  comparison  between  the  results  of  DABL 
and  TAPS  for  both  a  heated  and  unheated  axisymmetric  body.  The  coordinates,  potential 
flow  velocities  and  surface  temperature  distribution  of  the  ellipsoid  are  given  in  Appendix  B 
while  a  comparison  of  the  potential  flow  velocities  computed  by  DABL  and  TAPS  is  shown 
in  Figure  4. 

First,  for  the  unheated  case.  Figures  5,  6  and  7  compare  the  skin  friction  coefficient, 
boundary  layer  shape  factor  and  displacement  thickness  Reynolds  number  as  a  function  of 
arc  length  computed  by  DABL  and  TAPS  at  Rl  — 4.5  X  lO^/ft.  In  all  cases  there  is  good 
agreement  between  DABL  and  TAPS  with  slight  differences  occurring  near  the  ellipsoid  nose. 
Figures  8  and  9  compare  the  computed  values  of  the  disturbance  growth  rates  for  nondimen- 
sional  frequencies  cof  =  0.32  X  lO"*^  and  cuf  =  0.486  X  ICT*^.  respectively.  The  corresponding 
amplification  ratios  are  shown  in  Figure  10.  Both  the  growth  rates  and  amplification  ratios 
of  DABL  and  TAPS  are  in  close  agreement.  It  should  be  noted  that  the  higher  amplification 
ratios  calculated  by  DABL  are  consistent  with  the  results  given  by  vonKerczek  and  Groves 
(1978). 

For  the  heated  case,  the  surface  temperature  distribution  shown  in  Figure  1  1  was 
used  by  both  DABL  and  TAPS  for  the  ellipsoid  at  Rj^  =  4.49  X  1  O^/ft.  Figures  1 2,  1 3  and 
14  compare  the  skin  friction  coefficient,  shape  factor  and  displacement  thickness  Reynolds 
number  computed  by  DABL  and  TAPS.  The  computed  disturbance  growth  rates  and  ampli¬ 
fication  ratios  for  frequencies  cof  =  0.32  X  10“'^  and  cof  =  0. 18  X  10r4  shown  in 
Figures  15,16  and  1 7.  In  all  cases  DABL  produces  larger  growth  rates  and  amplification 
ratios  than  TAPS.  As  stated  by  vonKerczek  and  Groves,  linear  stability  theory  is  sensitive 
to  the  boundary  layer  shape  factor  and  Figure  13  shows  that  DABL  yields  a  slightly  larger 
shape  factor  than  TAPS  which  could  account  for  the  larger  growth  rates  and  amplification 
ratios.  Figure  17  also  compares  the  amplification  ratios  at  cof  =  0.22  X  10"“^  between  DABL 
and  TAPS  and  DABL’s  amplification  ratios  for  a  constant  surface  temperature  difference. 

AT  =  6.0C.  It  can  be  concluded  that  DABL  produces  more  conservative  results  as  compared 
to  TAPS. 
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Figure  4.  Potential  flow  velocity  distribution  of  a  9:1  ellipsoid. 
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Figure  6.  Boundary  layer  shape  factor  of  a  9:1  ellipsoid  -  Rl  =  4.5345455  X  10^/ft 
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Figure  7.  Displacement  thickness  Reynolds  number  of  a  9: 1  ellipsoid 
R,  =  4.5345455  X  lO^/ft. 
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Figure  8.  Computed  growth  rates  of  a  9: 1  ellipsoid  for 
Wf  =  0.486  X  10"^-  R,  =  4.5345455  X  loVt. 
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Figure  9.  Computed  growth  rates  of  a  9: 1  ellipsoid  for 
wf  =  0.32  X  1(H  -  Rl  =  4.5345455  X  lO^/ft. 
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Figure  10,  Comparison  of  the  computed  amplification  ratios  of  a  9: 1  ellipsoid 
for  Wf  =  0.486  X  1 0-4  and  Wf  =  0.32  X  1 0-4  -  =  4.5345455  X  1  O^/ft. 
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Figure  14.  Ksplacement  thickness  Reynolds  number  of  a  heated  9: 1  ellipsoid  - 
R  =  4.4901734  X  lO^t- 


GROWTH  RATE  a,  X  10 


0.20  0.30  0.40  O.50  0.60  0.70  0.80  0.90 

ARC  LENGTH,  S- ft 

Figure  15.  Computed  growth  rates  of  a  heated  9: 1  ellipsoid  for 
wf  =  0.32  X  1 0-'^  -  Rl  =  4.4901 1 34  X  lO^/ft. 
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Figure  1 6.  Computed  growth  rates  of  a  heated  9: 1  ellipsoid  for 
c>;f  =  0,18X  10-4- =  4  4901734  X  lO^/ft. 
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COMPUTER  PROGRAM  USAGE 
DEFINITION  OF  THE  INPUT  VARIABLES 


VARIABLE 

DESCRIPTION 

IBL 

[f  [BL  =  0,  calculate  the  boundary  layer  properties  and  perform  the  stability 
analysis 

[f  [BL  =  1.  calculate  the  boundary  layer  properties  to  separation  or  the  end 
of  the  body 

IBLTP 

tf  IBLTP  =  0,  do  not  read  or  write  boundary  layer  data  to  unit  8 

If  IBLTP  =  1,  write  boundary  layer  data  to  unit  8 

If  IBLTP  =  2,  read  boundary  layer  data  from  Unit  8 

IHEAT 

[f  [HEAT  =  0,  unheated  boundary  layer 
[f  [HEAT  =  1,  heated  boundary  layer 

I  POT 

[f  [POT  =  0,  calculate  potential  flow 

If  [POT  =  1,  calculate  potential  flow  and  punch  cards 

If  [POT  =  2,  input  potential  flow  according  to  Data  Table  4 

[START 

If  [START  =  0,  automatic  start  of  the  stability  analysis 

If  [START  =  I,  the  user  provides  information  -  Data  Table  2  -  for  starting/ 
restarting  the  stability  calculations 

ITWOD 

If  ITWOD  =  0.  axisymmetric  body 

If  ITWOD  =  1,  2-D  body’s  lower  surface 

If  ITWOD  =  2.  2-D  body’s  upper  surface 

[FIT 

[FIT  controls  the  amount  of  boundary  layer  velocity  and  temperature  profile 
data  used  in  the  stability  analysis  at  each  station 

If  I  FIT  =  1.  use  every  point  —  normal  option 

If  I  FIT  =  2,  use  every  second  point,  etc. 

[KFZ 

If  IKFZ  =  0,  stability  analysis  stops  when  the  disturbances  no  longer  experi¬ 
ence  amplification 

If  IKFZ  =  1,  stability  analysis  proceeds  regardless  of  amplification 

[V[SC 

If  I  Vise  =  1,  water  boundary  layer 

If  IVISC  =  2,  air  boundary  layer 

NP 

Number  of  body  offsets 

NP  <  160  axisymmetric  body 

NP^  101  2-D  body 

[BEG 

If  IBEG  =  0,  stability  analysis  autostart  begins  at  station  1 

If  IBEG  ¥=  0.  stability  analysis  autostart  begins  at  station  X(IBEG) 
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VARIABLE 


DESCRIPTION 


KF 


The  number  of  nondimensional  frequencies  for  which  the  amplification  ratio 
is  computed 


BEPS 


IfISTART  =  0,  3<KF<I0 
IfISTART=  1,  1  <KF<  10 

The  convergence  criteria  in  the  iteration  for  the  stagnation  point  in  the 
boundary  layer  code 

BEPS  =  0.0001  is  the  recommended  value 


BLDELS 

DELS 

TOL 

EPS 

X 

Y 

OMU 

SO 

RL 

ANORM 

R1 

U1 

JMAX 

ETA 

ALFRU 


Arc  length  increment  to  be  used  in  the  boundary  layer  calculation 

The  estimated  arc  length  spacing  at  which  the  instability  growth  rates  are  to 
be  computed.  The  program  varies  DELS  as  needed. 

The  accuracy  to  be  maintained  in  the  integration  of  the  stability  equation 

TOL  =  0.001  is  the  recommended  value 

The  convergence  criteria  for  the  eigenvalue  integration 

EPS  =  0.001  is  the  recommended  value 

Input  abscissa 

Input  ordinate 

Nondimensional  frequency  cof  for  which  the  amplification  ratios  are  to  be 
computed 

Arc  length  starting  values  for  the  stability  calculations  corresponding  to  OMU 

Reynolds  number  =  UgoL/r'oo 

Where  is  the  free  stream  velocity  (ft/sec); 

t'oo  is  the  kinematic  viscosity  (ft2/sec);  L  is  a  reference  length 

A  normalization  factor;  the  X,  Y  input  offsets  are  divided  by  ANORM 

Displacement  thickness  Reynolds  number  at  the  last  completed  station;  it  is 
input  only  when  restarting  the  stability  calculations  —  ISTART  =  1 

IfISTART  =  0;  R1  =0 

Potential  flow  surface  velocity  at  the  last  completed  station;  it  is  input  only 
when  restarting  the  stability  calculations  -  ISTART  =  1 

If  ISTART  =  0,  U1  =0 

The  number  of  points  across  the  boundary  layer  at  which  the  velocity 
profiles  are  computed 

JMAX  <201 

JMAX  =201  is  the  recommended  value 
The  scaled  edge  of  the  boundary  layer 
ETA  =  6.0  is  the  normal  input 

The  starting  estimate  or  the  value  at  the  last  completed  station  SO  for  the 
wave  number  (a^*);  it  cannot  be  zero  and  is  usually  less  than  0.5  for  a 
starting  estimate 


VARIABLE 

DESCRIPTION 

ALFIU 

The  starting  estimate  or  the  value  at  the  last  completed  station  SO  for  the 
disturbance  growth  rate  (a*);  a  value  of  zero  will  usually  suffice  as  a  starting 
estimate 

AMP 

The  value  of  the  amplification  ratio  at  the  last  completed  station  SO  or  zero 
for  a  starting  estimate 

ITT 

A  control  for  the  type  of  temperature  distribution 

If  ITT  =  0,  no  heating 

If  ITT  =  1 ,  a  step  function  temperature  distribution 

If  ITT  =  2,  a  linear  temperature  distribution 

If  ITT  =  3,  temperature  is  given  at  each  input  x  offset 

ITEMPl 

The  station  number  corresponding  to  the  offset  at  which  the  surface  tempera¬ 
ture  variation  is  to  begin 

If  ITT  =  3,  ITEMPl  =0 

ITEMP  2 

The  station  number  corresponding  to  the  offset  at  which  the  surface  tempera¬ 
ture  variation  is  to  end 

If  ITT  =  3,  ITEMP2  =  0 

TEMPI 

The  free  stream  temperature  in  degrees  Celsius 

DTEMP 

The  difference  in  temperature  AT  in  Celsius  between  the  free  stream  tempera¬ 
ture  and  the  maximum  body  temperature 

If  ITT  =  3,  DTEMP  =  0 

T 

The  body  temperature  corresponding  to  the  x  input  offsets 

UB 

The  potential  flow  surface  velocities 

XC 

The  scaled  abscissa  corresponding  to  the  maximum  2-D  body  diameter 

YC 

The  scaled  ordinate  corresponding  to  the  midpoint  of  the  maximum  diameter 

ANGO 

The  angle  of  attack  in  radians 

METHOD  OF  INPUTTING  DATA  TABLES 

The  input  data  tables  of  body  coordinates,  temperature  distribution  and  potential  flow 
are  input  as  shown  below. 


Body  Coordinates 


F  Format 


Card 

Columns 

1  -  10 

11-20 

1  I 

21-30 

iiiai 

31  -40 

41  -  50 

51-60 

Xl 

X2 

X3 

1 

•  •  • 

•  •  • 

•  «  « 

•  •  • 

.  •  . 

Xnp 

Yl 

^2 

Y3 

•  •  • 

•  •  • 

.  .  . 

•  •  • 

«  •  • 

^NP 
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Constant  Card  1 


Card 

Columns 

Format 

Variable 

1  -  20 

F 

RL 

21-30 

F 

ANORM 

31-40 

F 

R1 

41  -  50 

F 

U1 

Constant  Card  2 


Card 

Columns 

Format 

Variable 

1-3 

I 

JMAX 

11-20 

F 

ETA 

21-30 

F 

BEPS 

31-40 

F 

BLDELS 

41  -  50 

F 

DELS 

51  -  60 

F 

TOL 

61  -  70 

F 

EPS 

Data  Table  1 ;  Body  Coordinates  Y  =  f{x) 

The  maximum  number  of  entries  this  table  may  contain  is: 

a)  NP  =  160  for  an  axisymmetric  body 

b)  NP  =  101  for  a  two-dimensional  body 

Figure  18  illustrates  the  method  used  to  input  the  coordinates  for  an  axisymmetric 
body  while  Figure  1 9  should  be  used  for  two-dimensional  bodies. 


Figure  1 8.  Input  direction  for  the  axisymmetric  body  coordinates. 


Figure  19.  Input  direction  for  the  two-dimensional  body  coordinates. 


Data  Table  2:  Start/Restart  Information  for  the  Stability  Calculations 

If  ISTART  =  0,  do  not  input  this  table.  Tlie  number  of  entries  is  1  <  KF  <  1 0. 


Lara 

Columns 

Format 

Variable 

1  -  10 

E 

OMU(l)* 

11-20 

F 

SO  (I) 

21  -30 

F 

ALFRU(I) 

31-40 

F 

ALFIU(I) 

41  -  50 

F 

AMP(l) 

*Note;  Repeat  this  card  for  each  value  of  OMU  (I  =  1,KF). 

Constant  Card  3 

If  IHEAT  =  0,  do  not  input  this  card. 


Card 

Columns 

1 

Format 

Variable 

1 

I 

ITT 

3-5 

I 

ITEMPI 

7-9 

I 

ITEMP2 

11-20 

F 

TEMPI 

21  -  30 

F 

DTEMP 

Data  Table  3:  Body  Surface  Temperature  Distribution  T  =  f(x) 

If  IHEAT  =  0  or  ITT  ^  3,  do  not  input  this  table. 

The  number  of  entries  in  this  table  is  NP. 

Each  input  temperature  corresponds  to  each  input  x  offset  -  Data  Table  1. 

For  two-dimensional  bodies,  the  surface  temperature  distribution  follows  the  input 
scheme  of  Figure  1 9, 
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Data  Table  4:  Potential  Flow  Surface  Velocities  UB  =  f(x) 

The  number  of  entries  in  this  table  is  NP-1. 

The  potential  flow  velocities  are  at  the  midpoints  of  the  input  x  offsets  —  Data  Table  1 . 
For  two-dimensional  bodies,  the  potential  flow  velocities  follow  the  input  scheme  of 
Figure  19. 

Constant  Card  4 

If  IPOT  =  2  and  ITWOD  =  0,  do  not  input  this  card. 


Card 

Columns 

Format 

Variable 

1  -  10 

F 

XC 

11-20 

F 

YC 

21  -30 

F 

ANGO 

COMPUTER  PRINTOUT 


INCLUDED  ITEMS 

All  input  quantities  are  printed  out  to  allow  for  rapid  checking  for  errors.  The 
following  section  contains  the  definitions  of  those  output  variables  not  previously  defined. 
Below  is  a  list  of  the  items  included  in  the  program  output. 

1.  Control  parameters 

2.  Input  body  coordinates 

3.  Start/Restart  variables  for  the  stability  calculations  (ISTART  =  1) 

4.  Temperature  control  variables 

5.  Input  temperature  distribution  (IHEAT  =  1,  ITT  =  3) 

6.  Normalized  body  coordinates  and  temperature  distribution  (IHEAT  =  1) 

7.  Boundary  layer  parameters 

8.  Stability  autostart  variables  (ISTART  =  0) 

9.  Stability  calculations  (IBL  =  0) 

PRINTOUT  DEFINITIONS 


Boundary  Layer  Parameters 


VARIABLE 

S 


_ DEFINITION _ 

The  computed  arc  length  values  at  which  the  boundary  layer  calculations 
are  made 


UEF 

DUE 

RADIUS 


Potential  flow  surface  velocity 

Derivatives  of  the  potential  flow  surface  velocity  with  respect  to  arc 
length 

Body  radius 
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DEFINITION 


DELTN 


BETA  The  Hartree  parameter 

GAMMA  A  boundary  layer  scale  factor 

DELTN  Nondimensional  boundary  layer  thickness  times 

DI  Nondimensional  displacement  thickness  times  ^Rl 

H  The  boundary  layer  shape  factor 

CT  The  wall  skin  friction  coefficient  times  Vr^ 

CQ  Nondimensional  local  heat  transfer  coefficient  divided  by 

VISCR  The  ratio  of  fluid  viscosity  at  the  wall  to  that  of  the  ambient  fluid 

Note;  To  obtain  the  actual  value  of  the  nondimensional  boundary  layer  thickness,  displace¬ 
ment  thickness,  skin  friction  coefficient  and  heat  transfer  coefficient,  divide  DELTN, 
DI  and  CT  by  and  multiply  CQ  by  yjRj^- 


Stability  Automatic  Start 


DEFINITION 


ALPHAC 

OMEGAC 


Stability  Calculations 


The  arc  length 

The  critical  displacement  thickness  Reynolds  number 
The  displacement  thickness  Reynolds  number 
The  critical  value  of  a 

The  critical  frequency  associated  with  the  critical  displacement  thickness 
Reynolds  number 


DEFINITIONS 


The  arc  length 

The  potential  flow  surface  velocity 

The  boundary  shape  factor 

The  nondimensional  displacement  thickness  times 

The  wall  skin  friction  coefficient  times  ^R^ 

The  displacement  thickness  Reynolds  number 
A  boundary  layer  scale  factor 
The  body  radius 

The  nondimensional  frequency  —  cof 
The  real  part  of  the  complex  wave  number  -  a* 


VARIABLE 

DEFINITIONS 

ALFIU 

The  relative  growth  rate  of  the  disturbance  -  imaginary  part  of  a* 

AMP 

The  amplification  ratio  A/A^ 

ITER 

The  number  of  times  the  eigenvalue  problem  is  solved 

NSTEP 

The  number  of  integration  steps  of  the  Orr-Sommerfeld  integration 

Note:  The  DABL  code  solves  the  adjoint  Orr-Sommerfeld  problem  where  the  complex  con¬ 
jugate  of  a  is  the  eigenvalue.  Therefore  the  output  yields 

ALFRIU  =  Im  (a*)  <0  Damped  | 

=  0  Neutral  >  Disturbance 
>  0  Amplified  ) 
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SUMMARY 


Several  examples  illustrating  the  use  of  DABL  for  heated  and  unheated,  two- 
dimensional  and  axisymmetric  bodies  have  been  presented.  For  both  a  heated  and  unheated 
flat  plate,  DABL’s  automatic  starting  mode  calculated  nondimensional  frequencies  which  re¬ 
sulted  in  neutral  stability  points  falling  within  the  previously  defined  neutral  curves.  In 
general,  DABL  produces  unconservative  neutral  point  estimates  for  the  unheated  plate  and 
conservative  estimates  for  the  heated  plate.  However,  for  the  unheated  plate,  when  the 
Reynolds  number  and  disturbance  frequencies  are  reduced,  the  results  fall  closer  to  the 
neutral  curve.  Also,  the  amplification  ratios  calculated  by  DABL,  for  the  unheated  plate, 
compare  favorably  with  published  results. 

The  boundary  layer  properties,  growth  rates  and  amplification  ratios  computed  by 
DABL  for  an  unheated  9: 1  ellipsoid  compare  favorably  with  the  results  of  TAPS.  However, 
for  the  heated  9: 1  ellipsoid,  DABL  consistently  calculated  larger  growth  rates  and  amplifica¬ 
tion  ratios  than  TAPS  even  though  there  is  close  agreement  in  the  calculated  boundary  layer 
properties.  From  these  results  it  can  be  concluded  that  DABL  produces  conservative  results 
for  heated  axisymmetric  bodies. 

The  DABL  computer  code  can  be  considered  an  effective  tool  in  performing  para¬ 
metric  design  studies  for  both  heated  and  unheated  bodies.  For  unheated  bodies,  DABL  can 
be  used  for  predicting  the  location  of  boundary  layer  transition  by  the  e^  method.  With  its 
straight  forward  input  scheme,  DABL  is  easy  to  use  and  requires  a  minimum  of  user  interaction. 
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APPENDIX  A 


This  appendix  contains  an  abbreviated  computer  printout  of  the  automatic  start  - 
stability  analysis  of  a  heated  flat  plate. 
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APPENDIX  B 


This  appendix  contains  an  abbreviated  computer  printout  of  the 
stability  analysis  of  the  heated  9: 1  ellipsoid. 
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